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Abstract 

We examine the influence of the anisotropy of the superfluid energy gap and residual Fermi- 
liquid interactions in the triplet-correlated neutron liquid onto neutrino energy losses through 
neutral weak currents. The neutrino-pair emission caused by the pair breaking and formation 
processes and by the spin-wave decays is considered for the case of the 3 P2 pairing in the state 
with rrij = 0. The simple analytical formulae are obtained. A comparison with the previous results 
of the average-angle approach shows that the gap anisotropy leads to quenching of the neutrino 
emissivity caused by the pair recombination processes on about 15% and to substantial suppression 
of the spin-wave decays. Residual particle-hole interactions increase the energy losses in both the 
channels on about 5%. 
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I. INTRODUCTION 



Neutrino emission from a superfluid neutron liquid is currently thought to be the domi- 
nant cooling mechanism of the baryon matter, for some ranges of the temperature and/or 
matter density. One of the mechanisms leading the neutron star cooling consists on the 
recombination of thermally excited baryon BCS pairs into the condensate with emission of 
neutrino pairs via neutral weak currents It is generally accepted that, for temperatures 
near the associated superfluid critical temperatures T c , emission from pair breaking and for- 
mation (PBF) processes dominates the neutrino emissivities in many cases. Thisjdea has 
become the basis of the minimal cooling paradigm for superfluid neutron stars 

According to modern investigations 
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4|-|8(, the neutrino radiation caused by the pairing 
of baryons into singlet state is strongly suppressed, and the dominant neutrino emission, at 
the long-cooling epoch, occurs from the triplet-correlated condensate of superfluid neutrons 
which, as expected, exists in the superdense core of the star. Neutrino energy losses owing 
to the triplet PBF processes have been initially derived in Ref. ignoring the anomalous 



wea. 



£ interactions. The self-consistent approach to this problem was developed in Refs. 



10Ml2|. where the anomalous effective vertices are found as the solution of the Dyson's 



equations in the ladder approximation. Poles of the found analytical solutions [13 



indicate the existence of the undamped collective oscillations of the order parameter in the 
superfluid 3 P2 condensate of neutrons with rrij = 0. Some of the collective excitations with 
the excitation energy smaller than the superfluid energy gap are able to decay into neutrino 
pairs. This is, for example, the "so-called" normal-flapping mode which represents the 
nonunitary excitation with a nonzero average spin expectation value and can be termed by 
spin wave. The neutrino decay of such spin waves could be important for thermally-emitting 
neutron stars, presumably cooling through the combination of neutrino emission from the 
interior and photon cooling from the surface, the latter is responsible for their observed 



thermal emissions 15]. 



Exact solutions of the vertex equations are complicated by the anisotropy of the triplet 
order parameter. Therefore the preliminary analysis of the collective oscillations and the 
neutrino energy losses from the 3 P2 superfluid neutron liquid has been performed in the 



average-angle approximation replacing 



by its value averaged over angles 



13 



;he anisotropic energy gap in the quasiparticle energy 
, Q- The spin- wave energy, as found in this 
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approach, is proportional to the temperature-dependent energy gap amplitude A (T). 

However the more rigorous analysis with taking into account of the gap anisotropy has 
shown that the average-angle approximation is not justified for examination of the collective 
oscillations of the superfluid condensate. It was found that the gap anisotropy leads to 
a strong decreasing of the level energy along with a lowering of the temperature. One can 
expect that this is to suppress substantially the neutrino energy losses caused by the spin 
wave decays because the rate of neutrino losses is strongly dependent on the wave energy. 

In this paper we examine the effect of the anisotropy of the superfluid energy gap and the 
residual particle-hole interactions onto the neutrino energy losses from the 3 P2 superfluid 
neutron liquid with rrij = 0. We consider the neutrino emission caused by the PBF processes 
and by the spin- wave decays (SWD). 

The paper is organized as follows. The next section contains some preliminary notes. 
We discuss the effective strong interactions in the superdense neutron matter, the order 
parameter for the triplet pair-correlated system with strong interactions and introduce the 
notation used below. In Sec. Ill, we present the general expression for neutrino emissivity 
of the medium through neutral weak currents in terms of the imaginary part of the current- 
current correlator. In Sec. IV, we examine the anomalous vertices with taking into account 
of both the Fermi-liquid effects and the anisotropy of the superfluid energy gap for the case 
of 3 P2 pairing with rrij = 0, as is adopted in the minimal cooling paradigm jsj]. We also 
discuss the self-consistent response of the superfluid neutron liquid onto external neutrino 
field. Finally, in Sec. V, we evaluate the self-consistent neutrino energy losses from the PBF 
and SWD processes and compare them with the neutrino losses caused by modified urea 
processes and nn-bremsstrahlung. Section VI contains a short summary of our findings and 
the conclusion. A short Appendix contains some transformation of the previous results. 

In this work we use the standard model of weak interactions, the system of units h = c = 1, 
and the Boltzmann constant ks = 1- 

II. PRELIMINARY NOTES AND NOTATION 

The order parameter, D = D a p, arising owing to triplet pairing of quasiparticles in a 
degenerate Fermi system, represents a 2 x 2 symmetric matrix in spin space, (a,/3 =t, ■!•)• 
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Near the Fermi surface this matrix can be written as (see e.g. 18| ) 

D(n) = Ab(n)-&g , 



(1) 



where & = (o"i, (T 2 , 03) are Pauli spin matrices, and g = i&2- 

The spin-angle structure of the triplet condensate is defined by a vector b (n) in spin space 
which depends on the direction of the quasiparticle momentum p. The angular dependence 
of the order parameter is represented by the unit vector n = p/p which defines the polar 
angles (9, ip) on the Fermi surface, 

rii = sin 9 cos ip, n 2 = sin 9 sin <p, n-s = cos 9. (2) 

We assume that the (temperature-dependent) gap amplitude A is a complex constant (on 
the Fermi surface), and b (n) is a real vector which we normalize by the condition 

<6 2 (n)) = l. (3) 

Hereafter we use the angle brackets to denote angle averages, 
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The triplet pairing leads to the energy gap, Ab(n), in the quasiparticle spectrum which 
is in general anisotropic. We are mostly interested in the values of quasiparticle momenta 
p near the Fermi surface p ~ pf, where the quasiparticle energy is given by 



E p = ^jel + AW (n) , (5) 

with 

e p ~V F (p-p f ), (6) 

and Vf 1 is the Fermi velocity of the nonrelativistic neutrons. Here the fact is used that, 
in the absence of external fields, the gap amplitude A is real. 

The spin-orbit interaction among quasiparticles is known to dominate in the nucleon 



matter of a high density 19|, |20[ with the most attractive channel of interactions in the 3 P2 
state with s — l,j — 2,1 — 1 . In this case the order parameter in the superfluid system 
can be constructed with the aid of the set of mutually orthogonal complex vectors b m . (n) 



which generate standard spin-angle functions of the total angular momentum j = 2 and 
rrij = 0, ±1, ±2, so that 

b mj (n)<rg= ^2 (^af3\lm s j(llm s mi\2m j )Y ltmi (n) , (7) 
and are normalized by the condition 



These can be found in the form (See details in Ref. 14j ) 



b = y/T/2(-n 1 ,-n 2 ,2n 3 ) ,b x = -y^Jl{n^in Zl ni + in 2 ) 



b 2 = ^/3/4( ni + zn 2 ,z ni - n 2 ,0) ,b. mj = (-)"*' b^.. (9) 

In our approach it is necessary to distinguish the interactions in the channel of two 
quasiparticles from the interactions in the particle-hole channel. Since we are interested in 
values of quasiparticle momenta near the Fermi surface, p ~ Pf*i, the momentum transferred 
in the collision of two quasiparticles is of the order of 2pp. In this case the non-central spin- 
orbit and tensor interactions are most important in the superdense neutron matter. The 
attractive spin-orbit interaction which dominates in the channel of two quasiparticles is 



the basic reason of the neutron pairing 



20[. The most attractive channel corresponds to 



spin, orbital, and total angular momenta s — 1, I — 1, and j = 2, respectively, and pairs 
quasiparticles into the 3 P2 states with rrij = 0,±1,±2. The substantially smaller tensor 
interactions lift the strong paramagnetic degeneracy inherent in pure 3 P2 pairing and mix 
states of different magnetic quantum numbers. It is well known that the tensor interactions 
generate also a small 3 F 2 admixture to the ground state and markedly modify the 3 P2 energy 



gap 



21 



LH24J. 

The purpose of our study is however not the ground state but the linear response of 
the superffuid system onto external neutral weak currents. In this work we assume that 
the ground state of the superffuid neutron system and the magnitude of the energy gap at 
the Fermi surface are the known external parameters. Neutrino emissivity of the neutron 
system with a mixed 3 P2- 3 F2 superffuid condensate was examined in Ref. [12| . According 
to this work, incorporating of the small admixture of the 3 F 2 state requires of sophisticated 
calculations but does not affect noticeably the excitation spectra and the neutrino emissivity 



through neutral weak currents. Accordingly, throughout this paper, we neglect tensor forces. 



The pairing interaction, in the most attractive channel, can then be written as 



19] 



gT a ^ s (p, p') =V M) X>m»<^) Q/3 (g&b* m .(n'f) s , (10) 

where g = ppM* jis 2 is the density of states near the Fermi surface; and V(p,p') is the 
interaction amplitude. 

Consider now the interactions in the particle-hole channel. In our analysis, we shall use 
the fact that the Fermi-liquid interactions do not interfere with the pairing phenomenon 
if approximate hole-particle symmetry is maintained in the system, i.e., the Fermi-liquid 
interactions remain unchanged upon pairing. Near the Fermi surface, the Fermi-liquid effects 
are reduced to the standard particle-hole interactions: 

Q^ai^S (nil') = f (W) 5 a p5^ + g (nn') <x Q/3 <x 7( 5 . (11) 

For generality, Eq. (ITT)) should be supplemented with contributions from spin-orbit and 
tensor interactions. However, in uniform media, the momentum of particle-hole type excita- 
tions equals the transferred momentum q. We are interested in the medium response at the 
time-like momentum transfer, q < u ~ A. In this case the contribution from the noncentral 



interactions, which depend on the transferred momentum [241. |25I| . would be proportional to 
some power of q/pF "C 1 and vanish in the limit, q — y 0, which we are interested in (see 
below). 

III. NEUTRINO ENERGY LOSSES 

The emission of neutrino pairs is kinematically possible thanks to the existence of a 
superfluid energy gap, which admits the quasiparticle transitions with time-like momentum 
transfer k = (u,q), as required by the final neutrino pair: k — k\ + k 2 . We consider the 
standard model of weak interactions through neutral weak currents. After integration over 
the phase space of escaping neutrinos and antineutrinos the total energy which is emitted 
into neutrino pairs per unit volume and time is given by the following formula (see details, 



e.g., in Ref. |26|): 

_ G F J\f f*\ [«uQ{u-q) 



T^T duj I rf3 C(i)-i Im 10 (w ' q) {KK " " V) ' (12) 



5 



where 6 (x) is the Heaviside step-function; /x,u = 0,1,2,3 are Dirac indices; Af = 3 is the 
number of neutrino flavors; Gf is the Fermi coupling constant, and n^ ak is the retarded 
weak polarization tensor of the medium. 

In general, the weak polarization tensor of the medium is a sum of the vector- vector, axial- 
axial, and mixed terms. The mixed vector-axial polarization has to be an antisymmetric 
tensor, and its contraction in Eq. (j!2p with the symmetric tensor k^k u — k 2 g^ u vanishes. 
Thus only the pure- vector and pure- axial polarizations should be taken into account. We 
then obtain n(^ ak = CyU.y' + Cpjlpf, where Cy and Ca are vector and axial-vector weak 
coupling constants of a neutron, respectively. 

The Fermi velocity is small in the nonrelativistic system, Vp -C 1, and we can study the 
neutrino energy losses in the lowest order over this small parameter. We are interested in the 
time-like domain of the transferred energy and momentum, q < u, and u > A, in accordance 
with the total energy and momentum of escaping neutrino pairs. Since the transferred space 
momentum comes in the polarization functions of the medium in a combination qV F <C 
u), A, one can evaluate the polarization functions in the limit q = 0. Conservation of 
the vector current requires Uy U (u > 0, q = 0) = 0. This relation reflects the fact that 
the neutrino-pair emission through the vector channel of neutral weak currents is strongly 
suppressed in nonrelativistic systems [4]. Therefore we focus on the axial channel of the 
weak interactions, assuming H^ eak (u>, q = 0) ~ Cjll 1 ^ (lu, q = 0). Further simplification is 
possible due to the fact that, in the lowest (zero) order over the particle velocity, only the 
space component of the axial-vector vertex survives. This allows to write 

iC„ («, q = o) * cls^mi ( w , q = o) (is) 

with i, j = 1, 2, 3. 

The field interaction with a superfluid should be described with the aid of four effective 
three-point vertices. There are two ordinary effective vertices corresponding to creation of a 
particle and a hole by the field that differ by direction of fermion lines and - two anomalous 
vertices that correspond to creation of two particles or two holes. Accordingly, in graphical 
representation, the polarization tensor represents a superposition of loops incorporating the 
ordinary and anomalous vertices connected by the ordinary and anomalous Green functions, 
as depicted in Fig. [TJ We use the adopted graphical notation for the ordinary and anomalous 
propagators, G = — » — , G~(p) = — « — , = * « , and = * » . The analytic 
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FIG. 1. Graphs for the polarization tensor. Discarding the residual particle-hole interactions we 
show the ordinary vertices points. 

form of the quasiparticle propagators in the momentum representation can be written as 

G(p K ,p)=G(p K ,p)i, P^ =P(p K ,p)=F(p K ,p)bag, (14) 

= #t (-p K , -p) = gabF (-p K , -p) . (15) 
We define the scalar Green functions 

G(P.,P) = ^-^, F(p K ,p) = F(-p K ,-p) = -f^. (16) 

k p .t fx, p 

where p is the quasiparticle momentum, and p K = (2k + l)irT with k = 0, ±1, ±2, ... is 
the fermionic Matsubara frequency which depends on the temperature T. The quasiparticle 
energy is given by 

E p = yjel + A 2 6 2 (n) (17) 

with 

e p ~v F (p- p F ) . (18) 



IV. EFFECTIVE VERTICES AND POLARIZATION FUNCTIONS 

The anomalous effective vertices, which we denote as f W (n,w) and (n, tu), are given 
by infinite sums of the diagrams, taking into account the pairing interaction in the ladder 
app—onfl The ordinate vertices, r(n, u) , r (n, u) = r T (— n, u>), incor- 
)orating the particle-hole interactions can be evaluated in the random-phase approximation 



281 ] . This can be expressed by the set of Dyson equations symbolically depicted by graphs 
in Fig. 121 where the particle-hole interactions ffTT]) are shown by the shaded rectangle. Wavy 
lines represent the pairing interaction fllOp . The first diagram on the right-hand side of the 
first line is the three-point vertex of a free particle. In the nonrelativistic case, the bare 
axial- vector vertex is given by the spin matrices <x. 
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FIG. 2. Dyson equations for full ordinary and anomalous vertices. The particle-hole interaction is 
shown by the shaded rectangle. Wavy lines represent the pairing interaction. 

The analytic form of these equations, which are to be solved simultaneously with the gap 
equations, was repeatedly considered before. The equations can be substantially simplified 
due to the fact that we are interested in the processes near the Fermi surface. After a series 
of algebraic manipulations (see details in Ref. 11] ) the solution valid for q = can be 
obtained in the form: 

T« = ^B m >,)(<xb mj )s, (19) 



T< 2 > =^B m .(^) g(&b mj ) 



(20) 



and 



(21) 



where u v = 2nrj with rj = 0, ±1, ±2, ... is the bosonic Matsubara frequency, and (n, cu v ) = 
4> (— n, oj v ) and B mj (cu v ) are to satisfy the equations (we omit the uj v dependence of the 
functions): 



A Xr, 



1 f dn' ., 



J2 ^2oi (b mj xb) • B mi - 4PI„ 



with 



(22) 
(23) 



Xr, 



= ( (b^bm - b 2 ) A) - ii ( (b^b m ) Jo) - 2 ( (b^b) (bb mi ) Jo) . (24) 



The functions A (n) and Xq (uj v ; n) are given by 



A (n) = / de[^— tanh — - — tanh ] , (25) 



and 

XoK;n) = y^ ¥i ^ T -,tanh-, (26) 

with 



E= ^e 2 + A 2 b 2 (n) . (27) 
To solve the set of Eqs. f[2"2"j) . f[2"3l we expand <ft ( n ) i n spherical harmonics 

(n,w n ) = 4> hm (u n ) Y t , m (6, ip) , (28) 

Irn 

where 9 and ip are, respectively, the polar and azimuthal angles of the vector n, and write 

oo . 

s ( nn ') = E ^nk E ^ ( > *S. ^> > ( 29 ) 

1=0 » 

Inserting (J2SD and ([29]) into Eqs. ([22]) and ([23]) we obtain 

iu^ 1 



B mj = ^ xB ) Yi >™) ( 3 °) 



and 



>!.,,, — ^/~^K&lfl&mfl + oTjTT 7 ? ( TX~2 -V m ; F > m ' — Vl,m;l',m> j <?V,m'> (31) 



where we defined dimensionless quantities 



hm;V,m> = 4vr ^ — (*£ m 2tf (b m . X b) ) • (% m ,X i (b^ X b) ) , (32) 



and 

Vl,m;l',m> = 4vr (^* m PXolV,m') . (33) 

We will focus on the condensation into the state 3 P2 with rrij = which is conventionally 
considered as the preferable one in the bulk matter of neutron stars. In this case 

b (n) = b (n) , b 2 = l - (1 + 3n 2 3 ) (34) 

Since the function X (oj n , 6 2 ) is axial symmetric the integration over the azimuthal angle can 
be done in Eqs. f[32]) and f [33]) . By performing subsequent summation over rrij = 0, ±1, ±2 
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we find that only rrij = ±1 contribute, and X-i = Xi = X- 111 this way we obtain (the u v 
dependence is omitted): 

^l,m;l',m' = ~ ^mfiSm'flTj — ~ (35) 

* X 

with 

A K) = (b 2 (nl) X (co v , nl) P t (n 3 )> , (36) 
where Pi (x) are the Legendre polynomials, and 



Vl,m;l' ,m' &m,m'Vl,m;l' ,rrf 

In these notation we obtain the Eq. ( 13~T|) in the form 



>l,m 



v ^ 



3 ^ AA' 

2 4A 2 x 



+ T]l,m;l',m 4>V 



(37) 



(38) 



For this equation has only trivial solution (j)i >m ^o — 0. For <pi = 0/ jO Eq. ( 138]) takes 

the form 



with A, as given in Eq. (136]) and 



3 ^ AA' 

2 4A2 x 



+ ll,V <j>l> 



Iijl' M = ( p i ( n s) b 2 (n 2 3 ) x o (u v , nf) P v (n 3 )) . 



(39) 



(40) 



Since the <fi (n, oj) = 4> {— n , only even values of I contribute into the expansion fl28|) which 
takes the form 



i=even 



The Fermi-liquid parameters are not well known for a wide range of neutron densities 

, the parameters 
24j, at the Fermi 



we consider. Some known data allow to hope, however, that, in Eq. (1391) 
Qi/ (21 + 1) decrease rapidly for I > 2. For example, according to Ref. 
momentum pp = 1.7 fm _1 one has Qq = 0.842 while |$j2 — 0.043 8. Therefore we take the 
approximation g; = for I > 2 thus obtaining 

X 



(l + f o (Pl ))x + 2Qott 2 (Pl ) 
Inserting this expression into Eq. ( |22l) gives B = B-t 2 = 0, 



2 ' 



(42) 



B±i (u v ) = - 



iu v 1 3 



2A V 2 (1 + | 0O (m )) X + 2 0o ^ (b 2 IoY 



( A 



(43) 
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with 



X = (2 [-fijbtbx - (b*b) (bib)] Jo + (bjbi - P) A) . (44) 

In obtaining Eq. ( T4"3"j) we used the identities (3 = 7 0j o = (b 2 Io)- After the replacement 
iu n — >■ u + we obtain the analytic continuation to the retarded vertex for particles. The 
replacement iu n — > u> — zO gives the advanced vertex for holes. 
From Eqs. 1115} , (j20]) and (j2Tj) we find 



and 



f « (n, fi) = -/ (fi, [e* (o-bx) £ + e (d-b_ x ) $] , 
(n, O) = -f* (0, J/) [e*S (*b x ) + e£ (crb^)] 



with e = (1, i, 0), 



n 



2A (T) 

The functions / (Q, y) and <fi (fl, y) are given by 



w A(T) 

1/ ^ 



T 



ft(6 2 X ) 



and 



where 



2(l + | 0o (5 2 X o »x-2 0o fi 2 <& 2 X o )^ 



X 

(i + | fl o (VZo)) x - 2g ^ 2 <fe 2 x > 2 ' 



and 



dn 3 [6ft 2 (1 + nl) - 3nj (l - n 2 3 )] X (n 2 , Q, y) 

+ / rfn 3 (l-3n 2 )A(n 2 ,y)j , 
Jo 

i r 1 

(b 2 l ) = - J ^ 3 (l + 3r2 2 )X . 



(45) 
(46) 

(47) 
(48) 

(49) 



(50) 



(51) 



(52) 



(53) 



In obtaining Eq. (146]) we used the identity X (cu — iO) = Xq (uj + iO). 

Making use of the effective vertices ( 145|) - (14H|) and the propagators (|T4|) -(!T6l) one can 
calculate the axial polarization tensor by the diagrams of Fig. [TJ(See details in Refs. jl0l.ll3|). 
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For q = we obtain: 



nX (fi,y) 




(54) 



with b (n) = (-ni, -n 2) 2n 3 ) and 6 2 = § (1 + 3ri|), as given by Eq. (IMP . 

The imaginary part of IT^ that arises owing to the PBF processes originates from the 



part of the axial polarization tensor (1541 arises from the pole of the function / (fi, y) at 
Q = Q s < 6 min . This contribution describes the neutrino energy losses caused by the spin 
wave decays. The PBF and SWD processes operate in different kinematical domains, so 
that the imaginary part of the polarization tensor consists of two clearly distinguishable 
contributions, ImlT^ = ImIip Bf + ImIlg WD , which we now consider. 

V. NEUTRINO LOSSES 
A. PBF channel 

First we examine the PBF processes occuring at oj > 2A n = lb (n) A or, equivalently, at 
^ > 6 m i n . From Eq. (|54|) we obtain 



Inserting this expression into Eqs. (Fl3|) and (Fl2|) we calculate the contraction of Im IT; 
with symmetric tensor k^ky — k 2 g^ v to obtain 



function Xo(ri3]fl,y) at Q > b 



mm 



1 / y/2. A one more contribution into the imaginary 




(55) 





(56) 



where we use the local frame with Oz || b, and q\\ and q± are defined as 



2 1 / i~ N 2 2 2 2 

?|| = Jo (qb) ,q± = q ~ q\\ ■ 



(57) 
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The functions Xq, <fi and / are independent of the space momentum of the neutrino pair. 
Therefore the integral over d 3 q in Eq. ffl~2l) can be performed, and we obtain the neutrino 
emissivity in the form: 

y 7 n e 

'o 



32 



dft- 



Re, 



-ft Re/] <6 2 ImX ). (58) 



15tt 6 a - ■ J q e xp(2yft) - 1 

Imaginary part of the analytic continuation of the function fl26|) is given by 



Iml (tu + «0) = — sign (ft) tanh — 

4ft Vft 2 - & 2 V 2 



With 6 2 given in Eq. (Ell) one finds for ft 2 > 1/2: 



<& 2 ImX > 



tc 2 VQ 1 



(|f; a (i + 2tf) tmll ^)e(tf-i)e( 2 -tf) 



+ 



7T 



VE 1 
48~ft 



(2ft 2 + l) arcsin 



2ft 2 - 1 



- vW - 12 



x tanh My 1 (ft 2 - 2) 
The integration over dVL is now parted into two intervals: 



2/ft 



7 7 



Z^ 2 fl + 2ft 2 )ft 5 . 



ft 2 Re/ 



dft- 



ft E 



7^2 (e» n + 1) 



(2ft 2 + l) arcsin ■ 



Re< 



V2ft 2 - 1 

where we denote r = T/T c , and y = y (r); 

v/6 



ft 2 Im/ 



- V6ft 2 - 12 



1807T 4 



cIg 2 fPf m*mt^ 



2 er g 

A 



crrrs 



(59) 



(60) 



(61) 



(62) 



M and M* are the effective and bare nucleon masses, respectively; Tg c = T C /10 9 K. 

Equation (jSTJ) with / (ft, ?/) and (ft, j/) as given in Eqs. f )50|) and ([51]) improves previous 



results obtained in Ref. 



10], where the vertex function was evaluated in the BCS and 



average-angle approximation yielding 



(63) 
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Therefore before proceeding to the detailed analysis of the neutrino losses caused by the 
PBF processes, we examine the obtained equations for the above approximations. 
Replacing / (O) — > / av (fl) in Eq. ( 161]) leads to the following result 



e av (r) = € T 7 y 7 



2 f°° 
+ - / dfi 



'l + 2tt 2 )tt 5 



71 



1/V2 (e» n + 1) 

ft 5 



(e« n + 1)' 



(2f2 2 + l) arcsin 



>/3 



V2fi 2 - 1 



V6fi 2 - 12 



(64) 



The neutrino emissivity, as found in Ref. 



10( in the average-angle approximation, is written 



in the form of the two-fold integral over the Fermi surface and over the energy of escaping 
neutrino pairs. The integration over the Fermi surface in that expression can be performed 



o 



m 

Q- 
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FIG. 3. PBF neutrino emissivity in units eo versus the dimensionless temperature r = T/T c in 
various approximations marked in the plot. 

analytically. In the Appendix, we demonstrate that such the integration results in Eq. ( 164"|) 



14 



For numerical evaluation of the neutrino losses, with making use of Eqs. (I6T!) and (|64|) 
it is necessary to know the function y (r) which in general is to be found with the aid of 
the gap equation. For the case rrij = the function is well investigated. We can adjust for 
example the simple fit to y (r) = a/2 (1 — r) (0.7893 + 1.188/r), as suggested in Ref. js]. 

In Fig. [3] we show the neutrino energy losses computed according to the exact Eqs. ( JBTl 
l50l |5T|) in comparison with the energy losses in the average-angle approximation, as given 
by Eq. ( 164]) . The gap anisotropy leads to quenching of the PBF neutrino emissivity on 
about 15%, while the exchange particle-hole interactions with Qq = 0.842 slightly increase 
the neutrino energy losses. The neutrino emissivity which incorporates both the effects is 
less than the result obtained in the average-angle approximation on approximately 10 %. 



B. SWD channel 

The function X (n 3 ; fl,y) is real at Q < 6 min . From Eq. ( 155]) . ( 150]) and ( |5T]) we find in 
this domain 



(2a; 2 - 2q 2 - q\) b 2 l 

^ (l + | g o<5 2 Xo»x-2g fi 2 <6 2 Xo> 2 ' 

This function has a pole owing to existence of eigen oscillations of the condensate at the 
frequency u s satisfying the condition 

1 + ^ 0O (b 2 l )^j X ~ 2g n 2 s (b 2 l ) 2 = 0, (66) 

where all the functions are to be taken at Q = Q s , defined as 

«■ " 2KWY (67) 

As discussed in Ref. [ijj], the corresponding nonunitary oscillations look like the "normal- 
flapping" mode in 3 He-A jsoj]. 

Solution to Eq. (166]) can by found by assuming that the oscillation frequency of this wave 
is small, Q 2 1. This allows to neglect u 2 in the integrand of Eq. ( 126]) and write 

n t / \ A 2 f°° de E 
X (n 3 ; u 8 ) ~ I (n 3 ) = — J — tanh— . (68) 
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After this simplification the analytic solution to Eq. f|66|) can be written as: 



Jo dn 3 



3n|(l 



[l-3n 2 )A 



6 Jq 1 dn 3 (1 



Tin 



(69) 



where the function A is given by Eq. ( 1251) . and 



1 



b 2 I ) = 2 ^ dn 3 (l + 3nl)l . 



(70) 



In Ref. [U[ the right-hand side of Eq. ( l69l) was evaluated for g = in the average-angle 
approximation assuming that the anisotropic gap in the quasiparticle energy is replaced 
by its average-angle magnitude, b 2 A 2 (T) — > A 2 (T). Such approach results in Q 2 S — > 1/20 
or, equivalently, lu* v ~ A (T) / \/5 with a simple temperature dependence of the excitation 
energy only through the gap amplitude. However, the more accurate calculation has 
shown that, for the eigen modes, the average- angle approximation is valid only in the limit 
T — > T c . The gap anisotropy leads to a strong decreasing of the energy of the flapping mode 
at lowering of the temperature, and u s tends to zero when T — > 0, as shown in Fig. |H 




FIG. 4. The frequency uj s (r) of the spin wave at q = in units of 2A (r) versus reduced temperature 
r = T/T c . Solid curve is calculated for go = 0.842. Long-dashed curve corresponds to go = 0. 
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Near the pole Q ~ fl s we can approximate 



^weak {k^kv k 2 g^ u ) 



-2p {x - <fo 2 X » ( (2c 2 - 2gJ - g 2 ) b 2 l c 



(1 + fso (b 2 i )) dx/ow - 2 0o {m y (n-n s + io) 



(71) 



n=n a 



The displacement of the pole in a complex w-plane is chosen so that to obtain the retarded 
polarization. Such a displacement is equivalent to the presence of a delta-function imaginary 
part 



ImII weak faK - k 2 g^) 

2np (ffi 2 <6 2 X > - X ) ( ^ - 2q\ - q\) b 2 l £ 



(1 + | 0O <6 2 X » dx/dW - 20o (b 2 X o y 



■6(n-Q a ). 



(72) 



Q=Q S 



To evaluate this expression we expand the function Zq near the pole to obtain 



x ~ J + {n 2 - n 2 s ) h 



(73) 



with I , as given in Eq. ( 1681) . and 



~ A A f°° de E 
h = — / t^f tanh — -. 
4 ./_j; 5 2T 



(74) 



In this way we find 



d X /dn 2 ^^(2((l + nl)l c 



n\ (1 - n\) h 



(75) 



and 



6Q 2 ( (1 + nj) Jo ) - 3 (n 2 (1 - nl) J ) + < (l - 3n 2 ) A) 



The remaining calculations, similar to ones performed in the previous section, result in the 
following neutrino energy losses caused by the decay of spin waves: 



64 r T 7 y 7 Q 5 s 
eswD = £o— v6- 



37r exp (2yfi s 



l + |0o(6 2 Jo))5x/^ 2 -20 O (6 2 J c 



(76) 



3 \ " / J ~" \ 

In this expression the spin- wave relative frequency fl s (r) is defined by Eq. (1661) . The SWD 
neutrino emissivity, as obtained earlier in the average-angle approximation, can be obtained 
from this expression if to replace the relative frequency with a constant Q s — > 1/20. 
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It was already mentioned the gap anisotropy leads to a strong decreasing of the wave 
energy along with lowering of the temperature. This is to suppress substantially the neutrino 
energy losses caused by the spin wave decays because the rate of neutrino losses is strongly 
dependent on the wave energy. 
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FIG. 5. SWD neutrino emissivity from neutron matter versus the temperature T in logarithmic 
scale. The solid curve is calculated by Eq. (|6ip with taking into account the anisotropy of the 
energy gap and the residual particle-hole interactions in the spin-spin channel with the Landau 
parameter go = 0.842. Short-dashed curve is same but with go = 0. The long-dashed curve is 
calculated in the average-angle approximation and with go = 0, as given by Eq. (f64"J) . 

In Fig. O we show SWD neutrino emissivity from neutron matter versus the temperature 
T in logarithmic scale. We show also the neutrino emissivity as calculated in the anisotropic 
BCS approximation and with taking into account of both the anisotropy of the energy gap 
and residual particle-hole interactions. The Fermi-liquid effects lead to a minor (about 15%) 
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increase of the neutrino emissivity in the SWD channel, however this emissivity is small in 
comparison with the result obtained in the average-angle approximation 



C. Competitive neutrino processes 




FIG. 6. PBF and SWD neutrino emissivities versus temperature T in comparison with the modified 
urea and bremsstrahlung emissivities at kp = 1.7. 

The neutrino emissivity for the processes discussed above is plotted in Fig. El together 
with the modified urea and bremsstrahl ung emissivities with the suppression factors resulting 



from superfluidity as obtained in Ref. 31]. 



The emissivity from the PBF dominates everywhere below the critical temperature for the 
3 P2 superfluidity except the narrow temperature domain near the critical point, where the 
modified urea processes are more operative. Neutrino losses caused by the bremsstrahlung 
and SWD are less effective than in the PBF processes. 
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VI. SUMMARY AND CONCLUSION 



We have calculated the neutrino energy losses from the 3 P 2 ( m j = 0) superfluid neutron 
liquid with accurate taking into account of the anisotropy of the superfluid energy gap and 
minimal account of residual Fermi-liquid interactions. In our analysis we have used the 
integral expressions for anomalous three-point effective vertices derived in Refs. Q, ll^j ]. 
We have examined the neutrino energy losses through neutral weak currents caused by the 
pair breaking and formation processes and by the spin-wave decays. The corresponding 
neutrino emissivity in is given in Eqs. flBTj) and f!76|) . 



Earlier the neutrino losses have been calculated discarding the Fermi-liquid effects with 
using the average-angle approximation, where the anisotropic energy gap in equations for 
the anomalous vertices was replaced by its average-angle magnitude, b 2 A 2 (T) — > A 2 (T). A 
comparison of the results allows to make some inferences concerning validity of the average- 
angle approach. 

As shown above, the gap anisotropy leads to quenching of the PBF neutrino emissivity 
on about 15%, while the residual particle-hole interactions increase the PBF emissivity on 
about 5%. This difference is practically indistinguishable in the logarithmic scale in Fig. [61 
Up to accuracy about ten percents, one can neglect both the anisotropy of the effective weak 
vertices and the Fermi-liquid effects in the PBF processes and use the simple expression in 
the form of one-fold integral for practical estimates. 



However, the exact account of the anisotropy dramatically modifies the SWD neutrino 
losses. This fact has a simple explanation. As found in Ref. 10] , the relative spin- 
wave energy is constant in the average-angle approximation, £lf = l/Vb. Taking into 
account of the gap anisotropy leads to a strong temperature dependence of the energy of this 
collective mode [14J . The wave frequency diminishes and tends to zero at the lowering of the 
temperature, as shown in Fig. HI Since the rate of the neutrino losses is strongly dependent 
on the wave energy the anisotropy leads to a rapid decrease of the SWD energy losses along 
with lowering of the temperature. A comparison of different competitive processes presented 
in Fig. |6] allows to conclude that the dominant energy losses from the 3 P 2 superfluid neutron 
liquid are generated by the PBF processes except the narrow temperature domain near the 
critical point, where the modified urea processes are more operative. 
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Appendix A: Neutrino emissivity in the averege-angle approximation 



The neutrino emissivity, as obtained in Ref. 10] in the averege-angle approximation, is 
written in the form of the two-fold integral 



2 



/A poo 4 

t b2 l dx (TT^f' (A1) 



where z = \/ x 2 + b 2 y 2 . With the aid of the change x = y\/Q 2 — b 2 one can recast this 
expression to the form 

<- = ^GldM^M-Tyy- / £P [ *z ( j^f (A2) 

Further simplification is possible if to change the order of integration and write the emissivity 
in the form 

tav - -^GlClK PF M*TjrW dQ ° 5 / ^ e("-fe) (A3) 
15vr J J (l + ey Q yJ 4tt ^2 _ b 2 

where b 2 = (1 + 3n§) /2 and 9 (x) is the Heavyside's step-function. The integation over the 
Fermi surface gives 



+ ^[2Q 2 + 1) arcsin = - V6Q 2 - 12 j 9 (fl - V2j (A4) 



Inserting this expression into Eq. f ]A3j) we arrive to Eq. ( 164|) . 
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